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Abstract: We use AdS/CFT correspondence to study two-particle correlations in heavy ion colli- 
er . sions at strong coupling. Modeling the colliding heavy ions by shock waves on the gravity side, we 
^ , observe that at early times after the collision there are long-range rapidity correlations present in 
', the two-point functions for the glueball and the energy-momentum tensor operators. We estimate 
rapidity correlations at later times by assuming that the evolution of the system is governed by 
ideal Bjorken hydrodynamics, and find that glueball correlations in this state are suppressed at 
large rapidity intervals, suggesting that late-time medium dynamics can not "wash out" the long- 
^ ! range rapidity correlations that were formed at early times. These results may provide an insight 
on the nature of the "ridge" correlations observed in heavy ion collision experiments at RHIC and 
LHC, and in proton-proton collisions at LHC 
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1. Introduction 

In recent years it has been suggested that the medium of quarks and gluons produced in heavy 
ion collisions at RHIC goes through a strongly-coupled phase at least during some period of its 
evolution [1-4]. The Anti-de Sitter space/Conformal Field Theory (AdS/CFT) correspondence [5-7] 
is often used to study the dynamics of this strongly-coupled medium [8-29] : while it is valid only 
for J\f = 4 super- Yang-Mills (SYM) theory, there is a possibility that the qualitative (and some of 
the quantitative) results obtained from AdS / CFT correspondence may be applied to the real- world 
case of QCD. 

The main thrust of the efforts to study the djTiamics of the medium produced in heavy ion 
collisions using AdS / CFT correspondence has been directed toward understanding how (and when) 
the medium isotropizes and thermalizes [8-25,27,28]. The existing approaches can be divided into 
two categories: while some studies concentrated on the dynamics of the produced medium in the 
forward light-cone without analyzing the production mechanism for the medium [8-12,25,28], a 
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large amount of work has been concentrated on studying the colhsions by modehng the heavy ions 
with shock waves in AdSs and attempting to solve Einstein equations in the bulk for a collision of 
two AdSs shock waves [13-24,27]. Many of the existing calculations strive to obtain the expectation 
value of the energy-momentum tensor {T^y) of the produced medium in the boundary gauge theory 
[8-14,16,17,23-25,27,28], since this is the quantity most relevant for addressing the question of 
the isotropization of the medium. Other works address the general question of thermalization by 
noticing that it corresponds to creation of a black hole in the AdS bulk, and by constructing a 
physical proof of the black hole formation with the help of a trapped surface analysis [15, 18-22]. 

In this work we concentrate on a different observable characterizing heavy ion collisions: we 
study correlation functions in the produced expanding strongly-coupled medium. Correlation func- 
tions have become a powerful tool for the analysis of data coming out of heavy ion collisions, 
allowing for a quantitative measure of a wide range of phenomena, from Hanbury-Brown-Twiss 
(HBT) interferometry [30], to jet quenching [31] and Color Glass Condensate (CGC) [32]. In recent 
years a new puzzling phenomenon was discovered in the two-particle correlation functions measured 
in Au + Au collisions at Relativistic Heavy Ion Collider (RHIC) [33-35]: the experiments see cor- 
relations with a rather small azimuthal angle spread, but with a rather broad (up to several units) 
distribution in rapidity. This type of correlation is referred to as "the ridge". More recently the 
ridge correlations have been seen in high-multiplicity proton-proton collisions at the Large Hadron 
Collider (LHC) [36], as well as in the preliminary data on Ph + Ph collisions at LHC. 

Several theoretical explanations have been put forward to account for the ridge correlations. 
They can be sub-divided into two classes: perturbative and non-perturbative. Perturbative expla- 
nations, put forward in the CGC framework in [37-42], are based on the long-range rapidity corre- 
lations present in the initial state of a heavy ion collision due to CGC classical gluon fields [43-46] 
(see [47-49] for reviews of CGC physics). In [37,38] the authors invoke causality to argue that long- 
range rapidity correlation can only arise in the early times after the collision, since at later times 
the regions at different rapidities become causally disconnected. This is illustrated in Fig. |I|, where 
one can see that the gray-shaded causal pasts of two particles produced in the collision (labeled 
by arrows with momenta ki and /C2) overlap only at very early time (the red-shaded region). The 
authors of [38] then suggest that the late-time radial flow due to hydrodynamic evolution would lead 
to azimuthal correlations characteristic of the "ridge". Alternatively, the authors of [37,39-41] have 
identified a class of Feynman diagrams which generate azimuthal correlations in nucleus-nucleus 
collisions. 

The CGC correlations found in [37, 39-41] are based on purely perturbative small- coupling 
physics: however, it remains to be shown whether such perturbative dynamics contains large enough 
azimuthal correlations to account for all of the observed "ridge" phenomenon. In the scenario 
of [37, 38] CGC dynamics provides rapidity correlations, while azimuthal correlations are gener- 
ated by hydrodynamic evolution. As we have already mentioned, it is possible that the medium 
created at RHIC is strongly-coupled [1-4]: if so, hydrodynamic evolution would then be a non- 
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perturbative effect, malcing the scenario proposed in [37, 38] implicitly non-perturbative. Purely 
non-perturbative explanations of the "ridge" include parton cascade models [50], hadronic string 
models [51], and event-by-event hydrodynamic simulations [52]. The causality argument of [37,38] 
is valid in the non-perturbative case as well: one needs correlations in the initial state, either due to 
soft pomeron/hadronic strings interactions [50,51], or due to initial-state fluctuations [52], in order 
to obtain long-range rapidity correlations. In this work we will use AdS / CFT to address the theo- 
retical question whether long-range rapidity correlations are present in the non-perturbative picture 
of heavy ion collisions. At the same time we recognize that a complete understanding of whether 
the "ridge" correlations observed at RHIC and LHC are perturbative (CGC) or non-perturbative 
in nature is still an open problem left for future studies. 



r|=const 

T =const 

i 




Figure 1: Space-time picture of a heavy ion collision demonstrating how long-range rapidity correlations 
can be formed only in the initial stages of the collision, as originally pointed out in [37,38]. Gray shaded 
regions denote causal pasts of the two produced particles with four- momenta ki and k2 , with their overlap 
region highlighted in red. We have drawn the lines of constant proper time r and constant space-time 
rapidity r] to guide the eye and to underscore that late-time emission events for the two particles are likely 
to be causally disconnected. 

The goal of the present work is to study long-range rapidity correlations in heavy ion collisions in 
the strongly-coupled AdS/CFT framework. In order to test for the long range rapidity correlations 
observed in heavy ion collisions, we would hke to study the two-point function {tr F^j^{x) ti F^^{y)) 
of glueball operators tr (F^^^) right after the collision but before the thermalization. According to 
causality arguments of [37,38], one expects that the long range correlations in rapidity should occur 
at such early times. The choice of observable is mainly governed by calculational simplicity. The 
metric for the early times after the collision of two shock waves in AdSs was obtained in [14, 16, 17]: 
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after formulating the problem in Sec. ^ and presenting our general expectation for the answer in 
Sec. D, we use this metric to calculate the correlation function of two glueball operators in Sec. |[ 
(Since the glueball operator corresponds to the massless scalar field in the bulk, we compute the 
two-point function of the scalar field in the background of the colliding shock waves metric.) Our 
main result is that we do find long-range rapidity correlations in the strongly-coupled initial state, 
albeit with a rather peculiar rapidity dependence: the two-glueball correlation function scales as 

C{ki,k2) ~ cosh (4 Ay) (1.1) 

with the (large) rapidity interval Ay between them. We also show in Sec. ^ that the correlator 
of two energy- momentum tensors {T2 (x) T2 (y)) (with 1,2 transverse directions) exhibits the same 
long-range rapidity correlations. This should be contrasted with the CGC result, in which the 
correlations are at most flat in rapidity [37,38,40,41]. Indeed the growth of correlations with 
rapidity interval in Eq. ( |1.1| ) also contradicts experimental data [33-36] . Although we should not a 
priori expect an agreement between AdS/CFT calculations and experimental QCD data, we argue 
in Sec. ^ that inclusion of higher-order corrections in the AdS calculation along the lines of [17] 
should help to flatten out such growth, though it is a very difficult problem to demonstrate this 
explicitly. 

Using the causality argument of [37, 38] illustrated in Fig. |I] we also expect that after ther- 
malization the rapidity correlations should only be short-ranged. As a result, due to causality, 
the initial long ranged correlations can not be "washed away" and will be observed at later times. 
This explanation is analogous to the resolution of the 'horizon problem' in the cosmic microwave 
background radiation (CMB), where the observed near-homogeneity of the CMB suggests that the 
universe was extremely homogeneous at the time of the last scattering even over distance scales 
that could not have been in causal contact in the past. This problem was solved by assuming 
that the universe, when it was still young and extremely homogeneous, went through a very rapid 
period of expansion (inflation). As a consequence of inflation, different regions of the universe be- 
came causally disconnected, while preserving the initial homogeneity. The idea that we pursue here 
for the heavy ion collisions seems to be of similar nature. To verify the statement that late-time 
dynamics can not generate (or otherwise affect) long-range rapidity correlations we study glueball 
correlation again in Sec. ^ now using the metric found by Janik and Peschanski [8], which is dual 
to Bjorken hydrodynamics [53]. (This is done in the absence an analytic solution of the problem 
of colliding shock waves: despite some recent progress [14, 16, 17] the late-time metric is unknown 
at present.) Performing a perturbative estimate, we find that, indeed, only short-range rapidity 
correlations result from the gauge theory dynamics dual to the Janik and Peschanski metric. 

We summarize our results in Sec. |[ 
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2. Generalities and Problem Setup 



2.1 AdS/CFT Tools 

We start with a metric for a single shock wave moving along a light cone in the direction [8] in 
Fefferman-Graham coordinates [54]: 

ds^ = — {-2dx+dx' + ti{x')z^dx~^ + dxl + dz'^} (2.1) 



where 



= (2-2) 

Here x^ = x = (x^^x"^), dx\ = {dx^y + (dx^)^, z is the coordinate describing the 5th 

dimension such that the ultraviolet (UV) boundary of the AdS space is at 2; = 0, and L is the radius 

of the AdS space. According to holographic renormalization [55], (T {x~)) is the expectation value 

of the energy-momentum tensor for a single ultrarelativistic nucleus moving along the light-cone in 
the ^'''-direction in the gauge theory. We assume that the nucleus is made out of nucleons consisting 

of "valence gluons" each, such that (T (x^)) oc N^, and the metric (|2?l| ) has no A^-suppressed 

terms in it. The metric in Eq. (|2.1j ) is a solution of Einstein equations in AdSs: 

Ri.u + ^g^.u = o. (2.3) 

Imagine a collision of the shock wave ( p.l|) with another similar shock wave moving in the light 
cone x~ direction described by the metric 

ds^ = — {-2dx+dx' + t2{x^)z'^dx^^ + dx]_ + dz'^} (2.4) 



with 



t2{x^) ^ %{T,^Ax^)). (2.5) 

c 

Here we will consider the high-energy approximation, in which the shock waves' profiles are given 
by delta-functions, 

h{x-) = fii6{x-), hix^) = fi2S{x+). (2.6) 

The two scales /ii and /i2 can be expressed in terms of the physical parameters in the problem 
since we picture the shock waves as dual to the ultrarelativistic heavy ions in the boundary gauge 
theory [16,56] : 

^i~p+A?^f, /i2~P2A^<'. (2.7) 
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Here pf, ^^^e the large light-cone momenta per nucleon, Ai and A2 are atomic numbers, and Ai 
and A2 are the typical transverse momentum scales in the two nuclei [16]. Note that fj,i and /i2 are 
independent of Nc. 



The exact analytical solution of Einstein equations (|2.3|) starting with the superposition of the 



metrics (|2.1| ) and (|2.4| ) before the collision, and generating the resulting non-trivial metric after the 
collisions, is not known. Instead one constructs perturbative expansion of the solution of Einstein 
equations in powers of ti and t2, or, equivalently, /ii and /i2 [14,16,17,23,24]. At present the metric is 
known up to the fourth order in /x's [14,16,17], and also a resummation to all-orders in /i2 (/^i) while 
keeping /ii {^2) at the lowest order has been performed in [17]. The validity of the perturbatively 
obtained metric is limited to early proper times r = a/2x+ x~, see e.g. [17] (though indeed the 
fuUy-resummed series in powers of /ii, fi2 would be valid everywhere). Since here we are interested 
in the early-time correlations (and due to complexity of the /i2-resummed metric obtained in [17]), 
we limit ourselves to the 0{^i^2) nietric obtained in [14,16] in the Fefferman-Graham coordinates: 



[2 + G(i 



dx~^ dx + [ti{x )z'^ + F{x~^,x , z)] dx~ 



+ 



t2(x+)/ + F(x+,x",2) dx^'^ + [1 + H{x^,x~,z)] dx]_ + dz' 



The components of the metric at the order-/ii ^2 are 

F{x+,x-,z) = -\i{x+,x-)z'^ --d^_ho{x+,x~)z^ - —d^_hi{x+,x-)z^ 

6 16 

F(x~^ ,x~ , z) = —X2(x~^,x^)z'^ d'ihQ(x~^,x~) z^ dlhiix^ , x~) z^ 

6 16 



2 ho{x^, X") / - 2 /ii(x+, x')z^ + - ti{x') t2(x+) z^ 



G{x ,x , z) = 
H{x~^,x~,z) = hoi^x'^ ,x~) z"^ + hi{x^ ,x^) z^ 



where we defined [16] 

ho{x^ ,x^) 
Ai(a;+,x") 



d_ 

d. 



ti{x )t2{x^), hi{x^,x ) 



ti{x )t2{x'^) 



d 

ho{x^,x~), A2(x+,a;") = hQ{x^ ,x') 



9_ 



along with the definition of the causal integrations 



(2i 



(2.9) 



(2.10) 



rfx'+[...](x'+), ^[...](x-) 



dx'- [...]ix'-). 



(2.11) 
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Below we will calculate correlation functions of the glueball operators 

1 



J{x) 



(2.12) 



in the boundary gauge theory.^ According to the standard AdS/CFT prescription,^ the glueball 
operator is dual to the massless scalar (dilaton) field in the AdSs bulk [57] with the action 



q4> 



16 7r2L3 



d xdz \f^^ g dM(t){x, z) dN(f){x, z), 



(2.13) 



where M,N = {^,z), ^ = (0, 1,2,3) and correspond to 4D field theory coordinates, while z is 
the coordinate along the extra fifth (holographic) dimension. (As usual g = det gMN-) 
The equation of motion (EOM) for the scalar field is 
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0. 



(2.14) 



Using Eq. ( p.l4| ), the dilaton action evaluated on the classical solution can be cast in the following 
form convenient for the calculation of correlation functions: 



~ 16 7r2 L3 



d'^x [^/^g''''(f){x,z)^;,(f){x,z)] 



c 

16 ir^ 



2=0 



d^x (1>b{x) 



■dMx.z) 



2=0 

(2.15) 



In arriving at the expression on the right of Eq. ( |2.15|) we have used the metric in Eqs. (|2] 
( p.9|) , and ( p.lOj ), along with the standard assumption that the fields have the following boundary 
condition (BC) at the UV boundary, 0(x, z — )■ 0) = (psix), which allowed us to approximate near 
z = 



L 



10 



.10 



1 - -2«ti(x-)t2(a;^ 



10 



yW 



In arriving at Eq. ( p.l5P we have also demanded that^ 

■\/^ g^^ (f){x, z) dz(t>{x, z) — 7- as z — )■ oo. 



(2.16) 



(2.17) 



Define the retarded Green function of the glueball operator ( |2.12[ ) (averaged in the heavy ion 
collision background). 



Gr(Xi,Z2) 



-t9ixl-xl){[Jixi),J{x2)]). 



(2.18) 



^When defining the glueball operator we assume that in the boundary theory the gluon field A° is defined without 
absorbing the gauge coupling gyM in it, such that the field strength tensor F"^ contains the coupling gvM- 

^ Since A = 4, with A the confornial dimension of J{x) , the mass of the dual scalar field, ~ A (A — 4), is zero. 
•^As one can see later, our classical solutions satisfy this condition. 
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According to the AdS/CFT correspondence the contribution to the retarded Green function coming 
from the medium produced in the coUision is given by [58]^ 



= JT'lJf V (2.19) 

d(pB{Xi) d(pB{X2) 

where we subtract the action Sq of the scalar field in the empty AdSs space to remove the contribu- 
tion of the retarded Green function in the vacuum. The latter has nothing to do with the properties 
of the medium produced in the collision and has to be discarded. 

Later we will be interested in the Fourier transform of the retarded Green function 

GR{h,k2) = J d'xid^X2e-"''-'^'-"'^-^^GR{xuX2). (2.20) 
(We are working in the (— , +, +, +) metric in the boundary four dimensions.) 
2.2 Kinematics 

We have defined above = {k^ ± k^)/\/2, k = {k^, k'^), k± = \k\ and k"^ = k\ — 2k^ k~ = —m^. 
The particle rapidity, defined as, y = | In is a useful variable, since the rapidity difference 
between any pair of particles remains unchanged if we go from the center of mass frame to any 
other frame by performing a boost along the longitudinal direction, x^. On the other hand, when 
k^ ^ y ^ Up = lncot(6'/2), where yp is pseudorapidity, and 6 is the angle at which the particle 
emerges in the center of mass frame. Furthermore, defining = \/k\ + rn^ , we can rewrite the 
light-cone components of the momentum as: k'^ = m±e^/\/2 and k^ = m±e~y / \/2. In the case 
when kj_ ^ one has ^ kj_/2. 

Consider two identical on mass-shell particles with momenta ki = {ki,ki,k^) and ^2 = 
(A;J, k2 ,k2)- Assuming kf = k^ = —rn? and k^ = k2 = k, we obtain 

g2 = (k2 - kif = -2m^ -2ki- k2 = sinh^ ^ > , (2.21) 

where Ay = y2 — yi with yi and y2 the rapidities of the two particles. In case when k\ ^ m? and 
A|/ 3> 1, we have ^ 2k']_ cosh Ay ^ k\e^y . It is worth noting that the momentum difference is 
space-hke, since = > 0. 



''As was shown in [59,60] the right-hand side of Eq. ( ^.ig| ) contains contributions of both the retarded and advanced 
Green functions Gr and Ga- In the lowest-order calculation we are going to perforin here the Green functions are 



real, and, since ReG/j = YieGp = ReG^ (with Gp the Feynman Green function defined below in Eq. (2.34)), we do 



not need to address the question of disentangling the contributions of different wave functions to Eq. ( ^.1£ ) and will 



adopt the convention of [61,62] by calling the object in Eq. (2.19) a retarded Green function 
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2.3 Defining the observable in the boundary gauge theory 

Let us now specify the observable we want to calculate in the boundary gauge theory. Our primary 
goal is to study rapidity correlations using AdS/CFT. Ideally one would like to find correlations 
between produced particles. However, A/" = 4 SYM theory has no bound states, and, at strong 
coupling, it does not make sense to talk about individual supersymmetric particles. Therefore 
we will study correlators of operators, starting with the glueball operator defined in Eq. ( |2.12| ). 
One can think of the glueballs as external probes to A/" = 4 SYM theory (in the sense of being 
particles from some other theory in four dimensions), which couple to the gluons in A/" = 4 SYM, 
and therefore can be produced in the collision. Later on we will also consider correlators of the 
energy-momentum tensor T^,y, which should be also thought of as an operator coupling to a particle 
(in four dimensions) external to the A/" = 4 SYM theory. 

We start with the glueball production. To study two-particle correlations we need to find the 
two-particle multiplicity distribution 

(2.22) 



d?ki dyi d'^k2 dy2 



where ki, yi and A;^, y2 are the transverse momenta of the produced particles (glueballs) and their 
rapidities, and d'^k = dk^ dk'^. As usual we can decompose the two-particle multiplicity distribution 
into the uncorrelated and correlated pieces 

d'N d^N d^N d^N,,,, 



d^ki dyi d'^k2 dy2 d^ki dyi d'^k2 dy2 d^ki dyi d'^k2 dy2 



We are interested in computing the second (correlated) term on the right hand side of ( 2.23 ). We 
begin by writing it as 

^"'^ o,{mk,M)?) (2.24) 



d^ki dyi d'^k2 dy2 



where M{ki,k2) is the two-particle production amplitude. (Note that since we are primarily in- 
terested in rapidity dependence of correlators, we are not keeping track of prefactors and other 
coefficients not containing two-particle correlations.) 

For the correlated term in Eq. ( |2.23| ) the amplitude of inclusive two-glueball production in a 
heavy ion collision is 

M{ki,k2) oc y"dSic/^a;2e-*'=i-^^-*'=2-^^ (n|T{J(xi)J(a;2)} 1^1,^2), (2.25) 

which is a consequence of the LSZ reduction formula with T denoting time-ordering. Here \n) 
denotes an arbitrary state of the gauge theory which describes other particles which may be produced 
in a collision apart from the two glueballs. 
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The state l^i, A2) can be thought of as the vacuum in the presence of a source, with the source 
being the two nuclei with atomic numbers Ai and Consider first the expectation value of the 
energy-momentum operator {T^u) in a nuclear collision. According to the standard prescription we 
can write it as 

^ '-'^ ~ JVA,e'S[A]^^^^A]W4A] ^ 

where S[A] is the action of the gauge theory. For simplicity we only explicitly show the integrals 
over gauge fields in Eq. ( |2.26| ), implying the integrals over all other fields in the theory. The objects 
and are some functionals of the fields in the theory describing the two colliding nuclei. 

For instance, in the perturbative QCD approaches such as CGC, these operators are Wilson lines 
along x~ = and = light cone directions [43-46]. ^ 

Using operators and states in Heisenberg picture one can rewrite Eq. (|2.26|) as 

{T^uix)) = {AuA2\T^.{x)\AuA2). (2.27) 

Comparing Eq. (|2.27|) to Eq. ( p.26| ) clarifies the meaning of the |y4i,y42) state by demonstrating 
that the averaging in Eq. (|2.27| ) is over a state of vacuum in the presence of nuclear sources (which 
of course strongly disturb the vacuum). 

Using Eq. (|]25D in Eq. (|]2|) we obtain 



„, "^'^7:; , oc [ d''xid''x2d'x[d'x'2e-'''<^'-<^-'''<^'-^'^^ 
d^kidyid''k2dy2 J 

X ^(Ai,A2|T{J(x;) J(4)}|n) {n\T{J{x,)J{x2)}\A,,A2) (2.28) 

n 

where T denotes the inverse time-ordering and we have used the fact that J{x) is a hermitean 
operator. Summing over a complete set of states \n) yields 



OC J d'^xi d'^X2 d'^x[ d'^x'2 e 



d^ki dyi d'^k2 dy2 

X (Ai, A2I T { J(x;) J(4)} T { J(a;i) J{x2)] A2). (2.29) 

As one could have expected, in order to calculate two-particle production, we need to calculate 
a 4-point function given in Eq. ( p.29| ). This is, in general, a difficult task: instead we will use the 
following simplification. Begin by replacing the complete set of states \n) by states On{x) |y4i,y42) 
obtained by acting on our "vacuum" state 1^1,^2) by a complete orthonormal set of gauge theory 
operators (9„(a;), such that 



l">("l =1] j d*xOn{x)\AuA2) {Ai,A2\0i(x) (2.30) 



^Calculation of the expectation value of T^i, in CGC is reduced to perturbative evaluation/resummation of 



Eq. ( 2.26 ) (see e.g. [63] for an example of such calculation) 



- 10 - 



with the normahzation condition 

{A^,A2\OUy)On{x)\A^,A2) = 5nm5^''\x - y). (2.31) 



Using Eq. in Eq. (g]2|) we write 

d^kidyid^k2dy2 J 

X {A,,A2\T{J{x\)J{x'^)On{x)} |A,A2) (Ai,A2|r{Ot(x) J(Xi) J(X2)}|A,A2). (2.32) 

To evaluate Eq. ( p.32|) we have to insert all possible operators On{x) from the orthonormal set in 
it. Noting that J{x) is a gauge-invariant color-singlet operator, we conclude that only color-singlet 
On{x) would contribute. Also, since the final state in a scattering problem should be an observable, 
the operators C„ should be hermitean. The set of contributing C„(x)'s should therefore include 
the identity operator, J{x), Tfj_^{x), etc. 

As we will see below, since we are using the metric ( |2.8| ), which is a perturbative solution of 
Einstein equations to order /ii /i2, we can only calculate correlators to order /xi /i2 as well. Moreover, 
correlators which are independent of /ii and fi2 are vacuum correlators that we are not interested 
in. Correlators of order /ii or fi2 correspond to performing deep inelastic scattering (DIS) on a 
single shock wave similar to [61,62,64], and are thus not directly relevant to the problem of heavy 
ion collisions at hand. Thus in this paper we are only interested in correlators exactly at the order 
fii fi2 in the expansion in the two shock waves. Using such power counting it is easy to see that 
inserting the identity operator (normalized to one to satisfy Eq. ( |2.31| )) into Eq. (|2.32|) in place 
of OnS would give us a contribution of the order of fif which is the lowest order contribution 
to double glueball production. Inserting J{x) or T^iy(a;) into Eq. ( |2.32| ) instead of 0„'s would give 
zero. One can also see that replacing (9„.'s by higher (even) powers of J{x) or Tfj_y[x) (properly 
orthogonalized) in Eq. ( |2.32D would generate non-zero contributions, which are either higher order 
in 111 and /X2 or iV^-suppressed. We therefore insert the identity operator into Eq. (|2.32|) , which in 



the color space can be written as 1 = 5°"^ /Nc to satisfy normalization in Eq. ( p.31|) , and write 

oc / d^xi d^X2 d'x[ d% e-' ^i-(-i-i)-^ fc2-(x2-4) 



d^ki dyi d'^k2 dy2 

X -^{Ai,A2\T{J{x[)Jix',)}\Ai,A2) {Ai,A2\T{J{xi)Jix2)}\Ai,A2) [l + 0{l/N'^)] . (2.33) 

We have thus reduced the problem of two-glueball production to calculation of two-point correlation 
functions! Note that the prefactor of makes the counting right: since each connected 

correlator is order- A*"^, we see from Eq. (|2.33|) that the correlated two-particle multiplicity scales as 

as well, in agreement with perturbative calculations [37-41]. 

Defining Feynman Green function 

Gpiku k2) = J d^xi d^X2 e-^*^i-i-^*^2-x2 ^^^^ ^ | j^^^^ j^^^^l j^^^ ^2.34) 
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we can summarize Eq. (|2.33|) as 

(fN 1 
d^hdy,7k,dy, N!\GFih,h)f. (2.35) 

With the help of the retarded Green function 

Gnih, h) = j d^^i d^X2 e-*^i-i-^fc2-2 e{xl - x°) (Ai, A2I [J(xi), J{x2)] 1^1,^2) (2.36) 

and using the fact that at zero temperature [CfP = \Gr\^ [58], we rewrite Eq. ( p. 351) as 

d^N 1 

Therefore we need to calculate the two-point retarded Green function at the order fi2- This is 
exactly the kind of Green function one can calculate using the AdS / CFT techniques of Eqs. ( p. 191) 
and (^) . 

3. A Simple Physical Argument 

Before we present the full calculation of the two-particle correlations in AdS, we would like to give 
a simple heuristic argument of what one may expect from such a calculation. First of all, as we 
have noted already, we are going to expand the Green function, and, therefore, the bulk field (p 
into powers of /xi and /i2, stopping at the order-/ii /i2. To find the field at the order-/ii /i2 one 
has to solve Eq. (|2.14|) with the metric taken up to the order fj,i /i2. Since we are interested in the 
long-range rapidity correlations, our goal is to obtain the leading rapidity contribution from the 
calculation. Analyzing Eqs. ( |2.20| ), ( |2.15| ), and ( |2.19| ), one can conclude that the leading large- 
rapidity contribution comes from terms with the highest number of factors of light-cone momenta, 
i.e., from terms like k2 and /c|f A;^ (but clearly not from kf k^ = m']_/2 which is rapidity- 
independent). Taking M = = — in Eq. (|2.14 ) one obtains, among other terms, the following 



(leading-rapidity) contribution: 

(?(-^9!0o, (3.1) 

where 0o is the field at the order (/ii)° (yU2)° and g^2)^ is the metric at order-/ii fi2- Concentrating on 
order-2^ terms in the metric, which, according to holographic renormalization [55], are proportional 
to the energy-momentum tensor in the boundary theory, and remembering that the latter is rapidity- 
independent at order-/ii /i2 [14,16], we use energy-momentum conservation, df^T'^" = 0, which, in 
particular, implies that d- T + 9+ T~^~ = 0, to write 

9I2) = -^9(2)- (3-2) 
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Therefore Eq. ( |3.1| ) contains the term 



(3.3) 



which contributes to the field at order-/ii /i2, and, as follows from Eq. ( |2.20| ), resulting in a 
contribution to the retarded Green function in momentum space proportional to 



G 



R 



jX 9(2) (^2 



(3.4) 



with the Fourier transform of into momentum space. Since metric component at the 
order-/ii /i2 can not be rapidity-dependent [14,16], we see that Eq. (|3.4|) gives 



l|Aj/|»l 



~ e 



2(s/2-J/i) 



(3.5) 



Adding the ki -v^ /c2 term, arising from the g~^~^ component of the metric in Eq. (|2.14|) , we get 

^«l|A,|»i ~ cosh(2A|/). (3.6) 



Defining the correlation function 



U \ — rf^fcl dyi d?k2 dy2 
(Pkidyi cPk2dy2 



(3.7) 



and using Eqs. (|3.6|) and (|2.37|) to evaluate it we observe that at large rapidity intervals it scales as 



G{h,k2 



l|A3/|»l 



cosh (4 Ay). 



(3.8) 



Indeed the argument we have just presented relies on several assumptions: in particular it 
assumes that no other term in the metric would cancel correlations arising from the terms we have 
considered. To make sure that this is indeed the case we will now present the full calculation. The 
result of our simplistic argument given in Eq. ( |3.8[ ) would still turn out to be valid at the end of 
this calculation. 



4. Two-Point Correlation Function at Early Times 
4.1 Glueball correlator 

We now proceed to the calculation of the retarded Green function in the background of the metric 
(PD, following the AdS/CFT prescription outlined in Eqs. (F^ , (|2l5|) , and (^lOl) . 
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4.1.1 Bulk scalar field 



First we have to find the classical scalar field (p. Similar to the way the metric (|2.8|) was constructed 
in [16], we will build the scalar field cf) order-by-order in the powers of /ii and fi2, assuming /xi and 
/i2 are small perturbations. We would like to find the solution of Eq. ( |2.14D up to order 0{fj,ifi2)- 
For this we use the following expansion, 



(X, Z) = (j)o{x, Z) + Z) + z) + 02(a;, z) + . . . 



(4.1) 



where 0o ~ ^{1^2)1 ^a,b ~ C^(/Ui,2) and 02 ~ £'(^^1/^2)- We will use the standard method (see 
e.g. [61,62,64]) and demand that the boundary conditions at 2; — > are as follows: 



(j)o{x,Z^O) = (f)B{x), (j)a{x,Z^O) = (f)b{x,Z^O) = (f)2{x, Z ^ 0) 



0. 



(4.2) 



In this case the variation of the classical action with respect to boundary value of the field 0b 
required in Eq. (|2.19|) is straightforward. 

Using Eq. ( p.8|) in Eq. ( p.l4|) , and expanding the linear operator in the latter in powers of fii 
and /X2 up to order-/ii /i2 with the help of (p.9|) and ( 2.10 ), the EOM can be written explicitly in 
the form 



□5 + zUi dl + zU2dl + ^z^M 



{x,z) = 



(4.3) 



Taking into account that ti = ti{x ) and t2 = t2{x^), we give the following list of definitions: 



+ - 9, + □ 



4 , 



± 5 



1 



dx' 



(4.4) 



^4 



M = (D + z^] tit2Vi- ^Dtit2& 



Dtit2dl + 2 {D + 5z^) tit2 <9+a 



9_ 

+ 5 ti {d+ ^2) 5_ + 5 / 12 (a„ ti)d+ + 10 z^ ti t2d, + 2 z^ h t2 dl 



D = 96 



1 1 



+ 16 



1 1 

did! 



Substituting expansion (^TTI) into ([4.3|) , and grouping different powers of /ii and /i2 together we end 
up with the following set of equations, listed here along with their boundary conditions: 



□ 50O(X,Z)=O, 0o(x,Z ^ 0) = 0b(x) , 

□ 50a(x, Z) = -Z^ ti{x') dl 0o(x, Z) , (paix, Z -> 0) = 



□ 50fe(x, Z) = -Z^ t2{x^) dl 0o(x, Z) 



(x, 2; -> 0) = , 



(4.5a) 
(4.5b) 
(4.5c) 



□ 502(X,Z) = -Z^ti{x )dl4)b{x,z) - Z^t2{x^)dl4)a{x,z) - — M0o(x,z) , 02(x, 2: 0) = , 

(4.5d) 



12 
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where we also imply that all the solutions should be regular at z — )■ oo. 

To solve equations (|4.5|) it is convenient to introduce a Green function G{x, z, z') satisfying the 
equation 



U^G{x,z,z') = z'^6{z-z'). 

The Green function can be written as 

G{x,z,z') = z"" z'^ hiz^V^) K2{zy,/ni) , 

where 2;{<^>} = {min, max}{z, z'}. We can rewrite the inverse of Ds operator as 

1 dz' 

—f{x,z)= —G{x,z,z')f{x,z'). 
'-'5 Jo z 

Solving the first equation in ( [4 .51) we find 



1 



0o(a;,2) = -z'a^K2{z,/n-4)(j)B{x) 



From Eqs. ( |4.5b| ), ( [4.5c| ), and Eq. ( [4.5d| ) we have 



a{x,z) = [/ti9+0o] 



4>b{x,z 

1 



M 



x^z) = —zUidi — zH2di(j)o + —zU2di — zH^di(j)o ^ — 



□ 



□ 



12 



We have constructed the bulk scalar field which we need to find the correlation function. 
4.1.2 Glueball correlation function 



(4.6) 



(4.7) 



(4.8) 



(4.9) 

(4.10) 
(4.11) 



We can now calculate the retarded glueball correlation function using Eq. ( 4.11 ) in Eqs. ( 2.15|) , 
( p.l9| ), and (|2.2CI| ). It is straightforward to check that 



— dz G{x, z, z') 



J 2-S>0 



-z'^ U^K2{z Jul) 



(4.12) 



Using Eq. (^T^, along with Eqs. ( ^33]) , ( ^71^ , and ( CT) , we obtain 



(4.13) 



where 



F{k^,k2) =F,ik,,k2) + Fuik,,k2) 



(4.14) 
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with 



^2 ( Z 



/oo / , \ /"OO 

dz z^K2iz yjkl j / dz' z'^ K2 



(4.15) 



and 



Fuiki,k2) 



h.2 roo 



12 



^ dz z^K2 (^z 



96 



16 



~2 



1 

12 



l^l "-2 _|_ "'I "^2 



+ ^ A;JA;2" y rf^ z^ K2 \J~kl 



96 



96 



^A/-j^ Ay-j^ ^ Aj-j^ A/-|^ 

16^2 



+ 2" 



( ^A/fc2 



^Aj-j^ Ay-j^ y Aj-j^ Aj-j^ 



- — [2 k-^k:^ + fc+fcf + A;2-A;+] / dz z^K2[z \/ kf ) K2 { z sj k. 



+ ^ dz z^K2(^z^^ K,(^z^, 



We have defined 



(^-j^ — 2 /u"!^ H~ /ijj 



2 /u-j^ H~ /I'j 



(4.16) 



(4.17) 



with ki^_[_ = k2^± = k±. 

Before evaluating the obtained expressions further, let us comment on some of their features. 
First one may note that Eq. ( [4.13| ) contains a delta-function of transverse momenta of the two 
glueballs 5''^'*(fci+^2)- This demonstrates that at the lowest non-trivial order in fii and /i2 expansion 
(order-/ii (12) there will be nothing else produced in the shock wave collision apart from the two 
glueballs. Note that indeed a non-zero (T^;/) in the forward light-cone at the order- /ii /i2 found 
in [14, 16] indicates that a medium is created: however this strongly-coupled medium in the TV = 4 
SYM theory without bound states and confinement does not fragment into individual particles, 
and at late times simply results in a very low (and decreasing) energy density created in the 
collision, similar to the asymptotic future of Bjorken hydrodynamics dual found in [8]. Since in 
our calculation we have explicitly projected out two glueballs with fixed momenta in the final 
state, those two glueballs are all that is left carrying transverse momentum in the forward light- 
cone. (Leftovers of the original shock waves may also be present, though they would not carry any 
transverse momentum.) This picture is in agreement with the dominance of elastic processes in 
high energy scattering in the AdS/CFT framework suggested in [65]. 

Another important aspect of the result in Eqs. ( [4.15| ) and ([4.16| ) above is that the integrals over 
z and z' diverge for time-like momenta ki and k2, i.e., for kf = —m? and k\ = —m? corresponding 
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to production of physical glueballs of mass m. This result should be expected in A/" = 4 SYM 
theory: since there are no bound states in this theory, we conclude that there are no glueballs. 
Thinking of Bessel functions K2{z,^ 



fcf 2) Eqs. (|4.15|) and ([4.16|) as contributing to the wave 
functions of glueballs in AdSs space [66-68], we conclude that the lack of glueball bound states 
in the theory manifests itself through de-localization of these wave functions, resulting in "bound 
states" of infinite radii, both in the bulk and in the boundary theory (if we identify the holographic 
coordinate z with the inverse momentum scale on the UV boundary). Since the glueballs for us 
have always been some external probes of the TV = 4 SYM theory, we conclude that one has to 
define the probes by re-defining their wavefunctions. This can be accomplished, for instance, by 
introducing confinement in the theory, by using either the "hard- wall" or "soft- wall" models [66,69- 
77]. The inverse confinement scale would define the typical size of the bound states. Indeed such 
procedure would introduce a model-dependent uncertainty associated with mimicking confinement 
in AdS/CFT, but is unavoidable in order to define glueball probes. Besides, our main goal here is 
to calculate long-range rapidity correlations, which are not affected (apart from a prefactor) by the 
exact shape of the glueball AdSs wave functions. We therefore model confinement by modeling the 



glueball (external source) AdS wave functions by simply replacing K2{zy kl 2) K2{z A) in Eqs. 
( [4.15| ) and ( [4.16|) with A > related to confinement momentum scale. We then rewrite Eqs. ( [4.15|) 



and (|4. 16| ) as 



poo poo 

Fi{ki, k2) = dz z^ K2 (z A) / dz' z'^ K2 {z' A) 
Jo Jo 

X [{k^ktfh (QlZ<) K2 (QlZ>) + {ktk2 fl2 {Q2Z<) K2 {Q2Z, 



(4.18) 



p fOO 

Fn{k^M)~ / dzz'' [K2{zA)X 
J-^ Jo 



384 32 2' 



21 



4 2 

m I mi 



1 

12 



'^l "^2 _|_ "'I "^2 



-2n roo 



dz Z^ [K2{zK)f 



384 32 z^ 



384 32 z"' 



4 2 
m I mi 



+ 8z" 



1 r°° 

+ dzz' [K2izA)] 

c- rco 

- - [mi + k+k^ + k2kt] / dz z' [K2 {z A)y 

-L^ Jo 

+ / dzz^K2{zA) K,{zA) , 
-J Jo 



4 2 
m I mi 



(4.19) 



where we have also replaced all rapidity-independent factors with powers of either glueball mass m 
or 



m± = ^/kl + 



The contributions in Eqs. ( |4.18| ) and ([4.19| ) (or those in Eqs. ( |4.15|) and ( [4.16| )) to the retarded 
Green function ( p. 20 ) are shown diagrammatically in Fig. ^ in terms of Witten diagrams. There the 
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wiggly lines represent gravitons, while the dashed line denotes the scalar field. Crosses represent 
insertions of the boundary energy-momentum tensors of the two shock waves (/ii and (12)- Fi from 
Eq. ( |4.18| ) corresponds to the diagram on the left of Fig. 0, while Fu from Eq. ( |4.19| ) is given by 
the term on the right of Fig. 0. 




Figure 2: Diagrammatic representation of the correlation function calculated in this Section. 



It is important to note that the Green function given by Eqs. (|4.13|) , ( [4.14|) , (|4.18|) , and ( [4. 191) is 



indeed real, justifying the assumption we employed in stating that Eq. (|2.19|) provides us a retarded 
Green function. This can also be seen from the diagrams in Fig. |^ in which one can not cut the 
scalar propagator. The imaginary part of appears at higher order in ^1^21 when one has more 
graviton insertions in the scalar propagator, allowing for non-zero cuts of the latter. 

Let us now study the large-rapidity interval asymptotics of the obtained correlation function 
( [4.13| ). One can deduce from the kinematics described in Section |2.2| that 



ktk, = e-^y , k^kt = -^e^y , (4.20) 
such that when Ay = 1/2 — |/i 3> 1 we have 

Q\ = kl+ ml ^ ml e^^ , Ql = kl + 2ml e"'^^ ^ kl- (4.21) 
Therefore, the contribution from Eq. ( [4.18| ) becomes 

1*00 fOO 

Fi{kuk2)\^^^^^^ dzz'K2izA)J^ dz' z" K2 {z' A) {k^ k^fh {Q,z^) K2 (Qiz^) . (4.22) 
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To determine the large-Qi asymptotics of I2 {Qiz<) K2 (<5i-2>) note that, according to Eqs. 
and -2'^ h {Qiz<) K2 {Qiz>) satisfies 



z^z'^hiQiz^) K2{Qiz>) = z'H{z-z'). 



Hence, for Qi larger than the inverse of the typical variation in z we have 



z^z'^h{Qiz^) K^iQiz^ 
which, when used in Eq. ( [4 .221 ) yields 



v'3 



large Qi 



(4.23) 



(4.24) 



2048 (fci-fc+)2 512 ml 



I Aj/>1 



7 g^^io 7 A 



10 



(4.25) 



This result implies that the rapidity correlations coming from this term grow as e^^ at the early 
stages after the collision. 

On the other hand, the dominant contributions from the second term, Fii(ki, k2), are coming 
from the expressions in the second and the fourth lines of Eq. ( [4.19|) . They give 



Fii{ki, k2] 



I ~ ^ 

Ia|/>1 12 

5 

~ 12 



f^l 1^2 _|_ '^1 '^2 



dz / [K2izA)y 



{k+k- + k+k^) / dz z^ [K2 {zA)f 



384 32 z^ 



4 2 
m I mi 



(4.26) 



256 



21 Aio 



^0 

A2 42 A^ 
1-3 — + 

m 



^ 5 ml 



1280 m2 



21 Aio 



Combining Eqs. (|4.25|) and ( |4.26|) in Eqs. ( [1.13| ) and (|4.14|) we obtain 
64 A*"^ jji jj2 rn\ 



Gniki, k2 



21 



A 



10 



■6^'\k, + k,) 



,2 Ay 



A2 42 A^' 
mi 5 mi 



+ e 



Ay 



(4.27) 



which, dropping the second term in the curly brackets and using the + O — symmetry of the 
problem can be generalized to 



Gniki, k2 



l|A3/|»l 



— ^10 ^ (^i + ^2) cosh(2 Ay) 



21 



A^ 42 A^ 
mi 5 m I 



We thus conclude that at large rapidity separations 



GRiki,k2) 



|A3/|»1 



cosh (2 Ay) 



(4.28) 



(4.29) 
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in agreement with our estimate in Eq. ( p.6|) . 
Using Eq. ( p.37|) we conclude that 



~ cosh(4A?/) (4.30) 

|A?/|»1 



d'^ki dyi d?k2 dy2 

such that the two-glueball correlation function defined in Eq. ( p.7| ) scales as 



~ cosh(4A|/), (4.31) 

just like in Eq. ( p.8| ). We have demonstrated the presence of long-range rapidity correlations in case 
of strongly-coupled high-energy heavy ion collisions. The rapidity shape of the obtained correlations 
is very different from the "ridge" correlation observed experimentally at RHIC and at LHC [33-36]. 
It is possible that higher order in /ii and fi2 corrections would modify the rapidity shape of the 
correlation, putting it more in-line with experiments. We will return to this point in Sec. ^. 

Let us now pause to determine the parameter of our approximation. Until now we have, 
somewhat loosely, referred to our approximation as to an expansion in fii and /i2- However, these 
parameters have dimensions of mass cubed, and can not be expanded in. From Eq. ( [4.28|) we may 



suggest that the dimensionless expansion parameters are /^i/A^ and fi2/-^^, where A is the inverse 



glueball size. Thus our result in Eq. ( 4.28 ) dominates the correlation function only for 



(4.32) 

Since, as can be seen from Eq. (p.7p, /ii and /i2 are energy- dependent, these conditions limit the 
energy range of applicability of Eq. ( [4.28[ ). Eq. ( [4.32| ) also makes clear physical sense: since the 



metric ( |2.8| ) with the coefficients given by Eqs. ( |2.9| ) and Eq. ( 2.10 ) is valid only for early proper 
times T satisfying yUi_2 r"^ ^ 1 [16,17], we see that the glueballs have to be small enough, 1/A ^ r ~ 
/i^2^^5 to be able to resolve (and be sensitive to) the metric at such early times. 



Note also that the obtained Green function (|4.28| ) is not a monotonic function of m±: for 



m_L <^ A it grows with m± as m\, but, for m± ^ A it falls off as l/m^, peaking at = (28/5) A^. 
This translates into correlation function C{ki, ^2) first growing with m_|_ (and, therefore, k_i_) as 
for m± -C A, and then decreasing as l/m\ for m± ^ A. Similar non- monotonic behavior has been 
observed for "ridge" correlation experimentally [33-36] . While in CGC-based approaches [37-42] the 
maximum of the correlation function is given by the saturation scale Qs, and happens at k± ~ Qs, 
in our AdS/CFT case the maximum appears to be related to the inverse size of the produced bound 
state and its mass, such that it takes place at k± !^ a/A^ — m^. At this point it is not clear though 
whether such conclusion is a physical prediction or an artifact of the perturbative solution of the 
problem in the AdS space. 

In order to make a more detailed comparison with experiment one needs to improve on our 
AdS/CFT approach both by calculating higher-order corrections in /ii and /i2, and, possibly, by 
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implementing non-conformal QCD features, such as confinement, along the lines of the AdS/QCD 
models [66,69-77]. The latter modification would certainly change our glueball wave functions 
in the bulk, modifying the Bessel functions in Eqs. ( [4.18| ) and ( [4.19| ). However, while the use of 



AdS/QCD geometry may affect the m^-dependence of the correlation function ( 4.28|) , one may see 
from Eqs. ( [4.18| ) and ( ^.19| ) that such modification would not affect our main conclusion about the 
rapidity-dependence of the correlations shown in Eq. ( [4.31| ). The leading large-rapidity asymptotics 



of the correlation function (|4.31| ) results from the second term on the right-hand-side of Eq. (|4.19|) : 



modifying the glueball wave function would only change the coefficient in front of the rapidity- 
dependent part.^ Since the growth of correlations with rapidity does not reproduce experimental 
data [33-36], our conclusion is that the inclusion of higher-order corrections in ni and fi2 is the only 
possibility for AdS/CFT (or AdS/QCD) calculations to get in line with the data. 

4.2 Energy-momentum tensor correlator 

We have shown that there are long-range rapidity correlations in the glueball operator of Eq. ( p.l2|) 



in the strong-coupling heavy ion collisions. At the same time we would like to extend this statement 
to correlations of other operators. Energy-momentum tensor is a natural next candidate. Indeed the 
glueball operator (|2.12|) is a part of the energy- momentum tensor: hence correlations in {J{x) J{y)) 



probably imply correlations in (T^j,(x) T^y(|/)) as well. To show this is true we will present an 
argument below, largely following [78,79]. 

Consider a field theory whose dual holographic description is given by the metric of the general 
form 

ds"^ = g^^\ dx^ dx^ = /(x+, X-, z) dx\ + g^^,{x^, x", z) d^^ dC , (4.33) 

where x_ = {x^,x'^), dx\ = {dx^)'^ + ((ia;^)^,and C,'^ = {x^,x~,z). Now, consider small perturbations 
around the metric independent of x\a;^, qmn = gfiN + hMN{x^ ,x~ , z). We will work in the 
huz = gauge. The metric ( [4.33|) has a rotational 0(2) symmetry in the transverse plane. Under 
the transverse rotations one may naively expect {/in, hi2, h22} components to transform as tensors, 
{hoi, h^i, ho2, hs2} components to transform as vectors, and {hoo, hos, h^s} components to be scalars 
under rotations. However, rewriting the transverse part of the metric as 



hn hi2\ ^ f {hu + h22)/2 \ / (/^ii - /i22)/2 hu 

h2ih22j V {hn + h22)/2j \ /i2i -(/^ii - /i22)/2 



(4.34) 



we see that hu + h22 is also invariant under 0(2) transverse plane rotations. Hence the final 
classification of the metric components under 0(2) rotations is: {hn — h22, hi2} are in the tensor 
representation, {/loi, h^i, hQ2, ^2} are vectors, and {/loo, hos, /133, hn + /i22} are scalars [78, 79]. 

Using the above classification we see that we can assume that the only non- vanishing component 
of hMN is hi2 = h2i = hi2{x~^ , x~ , z) . It is in the tensor representation and, as can be seen 



^As the integrand in that term is positive-definite for any glueball wave function, the coeflieient ean not vanish. 
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with the help of Eq. ( 4.34 ), by rotating in the transverse plane we can always find a coordinate 
system in which hu — /122 = and hi2 = /i2i remains the only non-zero metric component in the 
tensor representation. Since all other components of the metric are in other representations of the 
0(2) symmetry group, they do not mix with hu in Einstein equations, and can be safely put to 
zero [78,79]. 

Substituting the metric qmn = S'mjv + ^A/iv( 

) with gf^j^ given by ( [4.33[ ) into Einstein 
equations (p.3|) , and expanding the result to linear order in hi2 we get [78,79] 

□ /ii2 - 2 d^hu + 2 ^-jf h,2 - y /^i2 = , (4.35) 

where 

□ = -L [V^g'''' dr,...] (4.36) 

and (5/)^ = g^^'^ OmI ■ Changing the variable from hi2 to h\ = hu/f, one can see that 
indeed satisfies the equation for a minimally coupled massless scalar [78,79]: 

ahl = 0. (4.37) 



Therefore, since our metric ( |2.8| ) falls into the category of Eq. ( |4.33| ), the analysis of Sec applies 



to the metric component hi. Defining the retarded Green function for the T2 components of the 
energy-momentum tensor (EMT) by 

G|*''^(fci,fc2) = J rf'xirf^X2e-^^^-^^-^^^-^^^(x?-x°)(Ai,A2| [T^{xi),T^{x2)] 1^1,^2) (4.38) 
we conclude that, similar to the glueball operator, 

Gr^(fci,fc2)||^^l>>, ~ cosh(2Ay). (4.39) 

Hence we have shown that the correlators of EMT operators exhibit the same long-range rapidity 
correlations as the glueball correlators. It is therefore very likely that such correlations are universal 
and are also present in correlators of other operators. 



5. Estimate of the Two-Point Correlation Function at Late Times 

Our conclusion about long-range rapidity correlations was derived using the metric ( |2.8| ) which is 
valid only at very early times after a shock wave collision. As discussed in the Introduction, we 
do not expect the interactions at later times to affect these correlations, since different-rapidity 
regions of the produced medium become causally disconnected at late times. To check that no 
long-range rapidity correlations can arise from the late-time dynamics one would have to calculate 
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the correlation function ( p.7| ) in the full metric produced in a shock wave collision including all 
powers of /xi and fi2- Since no such analytical solution exists, instead we will use the metric dual to 
Bjorken hydrodynamics [53] constructed in [8]. One has to be careful in interpreting the result we 
obtain in this Section: Bjorken hydrodynamics [53] is rapidity- independent, while there are reasons 
to believe that the medium produced in a shock wave collision would exhibit rapidity dependence, as 
indicated by perturbative solutions of Einstein equations done in [14,16,17]. Nonetheless, we expect 
that our calculation below would be a good initial estimate of the late-time rapidity correlations. 

The dual geometry corresponding to the perfect fluid was obtained by Janik and Peschanski 
in [8]. It can be written as 

/ \l/4 

where r = \/2x+x~ is proper time, 1] = ^ ln(x+/x~) is space-time rapidity, and -2/i(t) = ( ^) t^^^ 
(with Sq some dimensionful quantity) determines the position of the dynamical horizon in AdS^ 
such that the Hawking temperature is 

T(r)^^^^m"\-^':'. (5.2) 



TTZhir) TT V 3 

Unfortunately finding the glueball correlation function in Bjorken hydrodynamic state is equiv- 
alent to finding boundary-to-boundary scalar propagator in the background of the Janik-Peschanski 
metric ( |5.1| ), which is a daunting task: such propagator has not yet been found even for the static 
AdS Schwarzschild black hole metric. Instead, to estimate the correlations we will perform a per- 
turbative calculation. 

At late times, when t ^ Sq , assuming either that z is fixed or is bounded from the above 
(by let us say an infrared (IR) cutoff coming from the definition of the glueball wave function), we 
can consider the ratio u{t) = z/ Zhij) ^ 1 to be a small quantity. If so, we can expand the EOM 
for the scalar field (|2.14|) up to 0{u'^) obtaining 

□50(r, 7], x^, z) + [4 dl - 04] 0(r, r/, x^, 2;) = , (5.3) 
□50 = -z^d, + °40 , = \dr {jdA) - ^d^^ - Vi0 = (29+9_ - Vi) . 

Expanding the scalar field in the powers of u we write 

= 00 + 01 + • • • (5.4) 
where 0o ~ C (m°) and 0i ~ (9 (W^)- Substituting this back into Eq. (p.3|), we get 

C 4 

□5 00 = 0, 05 0, = -^A_ [492-04] 00. (5.5) 
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The solution for 0o was found above and is given in Eq. (O). We write the solution for 0i as 







3 r4/3 



3 r4/3 



□ 



4^0 



(5.6) 



where in the last step we neglected 9^, since a derivative like this generates 0(l/r^) corrections (at 
fixed u), which were neglected in constructing the original metric ( p.l| ) and are thus outside of the 
precision of our approximation. We are now ready to calculate the retarded Green function. Using 
Eq. in Eqs. ( pJ9D , ( pJfSD , and (^[20]), and employing Eq. ( ^:i2D yields 



lo{i/4) 



iV,^gom^ 
24 



/■oo 

^^(^1 + ^2) / dzz^K2 
Jo 



zJkf K2 zJk. 



X 



1 1 

c/x'*' dx" e'"""^ {k-+k-)+ix- ik++k+) 

7-4/3 



(5.7) 



where we have replaced kf and k^ with — everywhere except for the arguments of the Bessel 
functions. The integrals over and x~ in Eq. ( ^.71) run from to cxd since the matter only exists in 



the forward light-cone. (On top of that the metric ( |5.1|) is valid at late times only, for -u ^ 1, such 
that the actual and x~ integration region should be even more restricted, possibly suppressing 
the correlations we are about to obtain even more.) 



Just like in the case of the early times considered in Sec. ^]T|, the integral over z in Eq. (|5.7| ) 

we recognize 



is divergent for time-like momenta ki and ^2. Similar to what we did in Sec. 4.1 



the Bessel functions in Eq. (|5.7|) as the glueball wave functions in the bulk, which need to be 
modified to refiect the finite size of glueballs, which do not exist in Af = 4 SYM theory. Replacing 
K2{zJk'l2) K2{z K) in Eq. (|5?7|) and integrating over z yields 



Gl'{k^,k2) 



15 A6 

Evaluating the integrals left in Eq. (|5.8| ), 

dx'^ dx~ e*^"*" (^r+'=2')+«^" i^t+^t) 



+^ 1 



r 



4/3 ■ 



(5.^ 



TV 



(2 x+ x-)2/3 (kt + k^y/^ik^ + A;2-) 1/3 



(5.9) 



where 



N 



r2 (1) e^-/3 

227^ 



we obtain 



G'r^(/Ci,/C2)|o(i/^4) 



N 



15 A6 



rn'l^ (1 + cosh Ay) 1/3 



(5.10) 



(5.11) 
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The corresponding two-glueball correlation function scales as 



C^Uki,k2)\,^ 7T. • (5.12) 

^ ''^\^y\»^ mf (coshAy)2/3 ^ ^ 

We conclude that rapidity correlations coming from the AdS dual of Bjorken hydrodynamics are 
suppressed at large rapidity interval, at least in the perturbative estimate we have performed. This 
result appears to agree with the causality argument [37,38] making appearance of long-range rapidity 
correlations unlikely at late times. Moreover, the locality of C^^ in rapidity suggests that late-time 
dynamics is not likely to affect long-range rapidity correlations coming from the early stages of the 
collision: hydrodynamic evolution can not "wash out" such long-range rapidity correlations. 

Note that the complete momentum space two-glueball correlation function receives contribu- 
tions from all regions of coordinate space, i.e., from all xi and X2- In Sec. § we have calculated 
the contribution arising from early proper times, while here we have estimated the late-time con- 
tribution. One may expect that in the complete result the two contributions coming from different 
integration regions would simply add together: in such case clearly the early-time contribution in 
Eq. ( |4.31|) would dominate for large rapidity intervals, leading to long-range rapidity correlations 



arising in the collision. 



6. Summary 

Let us summarize by first restating that we have found long-range rapidity correlations in the initial 
stages of strongly-coupled heavy ion collisions as described by AdS/CFT correspondence. We expect 
that due to causality the correlations would survive the late-time evolution of the produced medium, 
though one needs to have a full solution of the shock wave collision problem to be able to verify this 
assertion. The long-range rapidity correlations may be relevant for the description of the "ridge" 
correlation observed in heavy ion and proton-proton collisions [33-36]. Indeed "ridge" correlation 
is characterized not only by the long-range rapidity correlation, but also by a narrow zero-angle 
azimuthal correlation between the triggered and associated particles. As was suggested in [37, 38] 
such azimuthal correlation may be due to the radial flow of the produced medium. The advantage 
of the AdS/CFT approach to the problem is that the full solution to the problem for a collision 
of two shock waves with some non-trivial transverse profiles would have radial flow included in the 
evolution of the dual metric, and would be able to demonstrate whether radial flow is sufficient to 
lead to the "ridge" phenomenon. Indeed such calculation appears to be prohibitively complicated 
to do analytically at the moment. 

The correlations we found grow very fast with rapidity interval, as one can see from Eq. (|3.8|) , 
while the experimentally observed correlation [33-36] is at most fiat in rapidity. This result may 
lead to the conclusion that the initial stages of heavy ion collisions can not be strongly-coupled, 
since this contradicts existing observations. At the same time, it may happen that higher-order 
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corrections in /ii and /i2 would affect this rapidity dependence, ffattening the resulting distribution. 
On yet another hand, such higher-order corrections become important at later times, and eventually 
causality may prohibit further late-time modification to the long-range rapidity correlations. More 
work is needed to clarify this important question about the rapidity- shape of the correlations coming 
from the solution of the full problem in AdS. 

Assuming that the issue of rapidity shape would be resolved, we would also like to point out that 
/ct- dependence of obtained correlator ( [4.28|) closely resembles that reported in the data [33-36] : it 
starts out growing with kr at low-fcy, and, at higher fcy, it falls off with kr- The location of 
the maximum of the correlator in our case was determined by the mass and size of the produced 
particles, and was thus energy-independent. It is possible that the solution of the full problem, 
resumming all powers of fii and fi2 would lead to the maximum of the correlation function given 

1 /3 

by fii 2, which in turn would be inversely proportional to the thermalization time [14,22], thus 
providing an independent way of measuring this quantity. Again more research is needed to explore 
this possibility. 
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